home *** CD-ROM | disk | FTP | other *** search
/ Amiga Format CD 46 / Amiga Format CD46 (1999-10-20)(Future Publishing)(GB)[!][issue 1999-12].iso / -in_the_mag- / reader_requests / scilab / man / man-part1 / cat7 / yulewalk.7 < prev    next >
Text File  |  1999-09-16  |  1KB  |  67 lines
  1.  
  2.  
  3.  
  4. yulewalk(1)                    Scilab Function                    yulewalk(1)
  5.  
  6.  
  7.  
  8.  
  9.  
  10.  
  11. NAME
  12.   yulewalk -  least-square filter design
  13.  
  14. CALLING SEQUENCE
  15.   Hz = yulewalk(N,frq,mag)
  16.  
  17. PARAMETERS
  18.  
  19.   N         : integer (order of derired filter)
  20.  
  21.   frq       : real row vector (non-decreasing order), frequencies.
  22.  
  23.   mag       : non negative real row vector (same size as frq), desired magni-
  24.             tudes.
  25.  
  26.   Hz        : filter B(z)/A(z)
  27.  
  28. DESCRIPTION
  29.   Hz = yulewalk(N,frq,mag) finds the N-th order iir filter
  30.                     n-1         n-2
  31.         B(z)   b(1)z     + b(2)z    + .... + b(n)
  32.   H(z)= ---- = ---------------------------------
  33.                   n-1       n-2
  34.         A(z)    z   + a(1)z    + .... + a(n)
  35.   which matches the magnitude frequency response given by vectors frq and
  36.   mag.  Vectors frq and mag specify the frequency and magnitude of the
  37.   desired frequency response. The frequencies in frq must be between 0.0 and
  38.   1.0, with 1.0 corresponding to half the sample rate. They must be in
  39.   increasing order and start with 0.0 and end with 1.0.
  40.  
  41. EXAMPLE
  42.   // f=[0,0.4,0.4,0.6,0.6,1];H=[0,0,1,1,0,0];Hz=yulewalk(8,f,H);
  43.   //fs=1000;fhz = f*fs/2; //xbasc(0);xset('window',0);plot2d(fhz',H');
  44.   //xtitle('Desired Frequency Response (Magnitude)')
  45.   //[frq,repf]=repfreq(Hz,0:0.001:0.5);
  46.   //xbasc(1);xset('window',1);plot2d(fs*frq',abs(repf')); //xtitle('Obtained
  47.   Frequency Response (Magnitude)')
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  
  66.  
  67.